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THE CALCULATION OF STEADY NON-LINEAR TRANSONIC FLOW 


OVER FINITE WINGS WITH LINEAR THEORY AERODYNAMICS 

By Atlee M. Cunningham, Jr. 

Fort Worth Division of General D 3 mamics 

SUMMARY 


This report docxaments a study to establish the feasibility 
of calculating steady mean flow solutions for non-linear tran- 
sonic flow over finite wings with a linear theory aerodynamic 
computer program. The methodology is based on independent solu- 
tions for upper and lower surface pressures that are coupled 
through the external flow fields. Two approaches for coupling 
the solutions are investigated which include the diaphragm and 
the edge singularity method. The final method is a combination 
of both where a line source along the wing leading edge is used 
to account for blunt nose airfoil effects; and the upper and 
lower surface flow fields are coupled through a diaphragm in 
the plane of the wing. An iterative solution is used to arrive 
at the non-uniform flow solution for both non-lifting and lift- 
ing cases. Final results for a swept tapered wing in subcritical 
flow, show that the method converges in three iterations and 
gives excellent agreement with experiment at oC =0° and 2°. 
Recommendations are made for development of a procedure for 
routine application. 


SYMBOLS 


[A] 

C-} 


p 


aerodynamic matrix 
pressure series coefficients 

pressure coefficient 

local wing chord, meters 
frees tream and local Mach niimber 
pressure. Newtons /me ter ^ 
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2 

d5mamic pressure. Newtons /me ter 

leading edge "radius”, meters 

frees tream velocity, meters /second 

velocity components in the x,y,z directions, 
meters /second 

total velocity, meters /second 

normal wash on surface Zs 

Cartesian coordinates, meters 

surface thickness coordinate, meters 

coordinates of the surface induced by 
the leading edge line source, meters 

upper and lower surface coordinates, meters 

meters 

meters 

angle of attack 
angle of surface Zs 

angle of the total velocity component, V 



INTRODUCTION 


Due to the strong coupling between the oscillatory and 
steady flow fields in the transonic speed range, it is no longer 
possible to calculate unsteady aerodynamic forces that are inde- 
pendent of the mean flow field characteristics. A completely 
adequate transonic aerodynamic method should account for the 
effects of surface thickness, twist, and camber distribution, 
mean angle of attack, wing-body and wing -wing interference, 
mixed flow with imbedded shocks, and boundary layer and boundary 
layer-shock interaction. Since these effects are embodied in 
the mean flow fields, they can be accounted for in the unsteady 
forces in a linear fashion as long as the oscillatory solutions 
are obtained for small perturbations about the mean flow. Thus, 
the nonlinearity can be contained within the steady solution 
which permits linearization of the unsteady solution (ref. 1). 

The coupling between steady and unsteady flow fields in the 
transonic regime has therefore made the calculation of steady 
state solutions an integral part of the unsteady solution pro- 
cesses used for transonic flutter and d}mamic response analyses. 
As a result, the usual Mach-altitude envelope for flutter clear- 
ance must be expanded to a three dimensional Mach-altitude-angle 
of attack envelope. Mean flow field data could be obtained from 
either experiment or finite difference solutions. From a prac- 
tical applications standpoint, however, it would be desirable 
to have the ability to calculate this data with the same program 
which is used to calculate the oscillatory solutions. Such a 
computer program was developed under Contract NASl-12399 (ref. 1) 
in which the piecewise linearization concept was applied to 
steady flow as well as unsteady flow. It appears that this 
computer program has the potential for developing an economical 
steady solution which will be suitable as a source of flow field 
information for calculating oscillatory solutions. 

This report docxaments a study on some candidate modifica- 
tions to the transonic computer program for providing the 
capability to calculate steady mean flow solutions. The 
methodology is based on independent solutions for upper and 
lower surface pressures which are coupled through the external 
flow fields. The two methods investigated are the diaphragm 
and the edge singularity methods. The final method is a com- 
bination of both where a line source along the wing leading 

edge is used to ac c ouiit ~f of ' hluh t no se” ai t f o i t^e f f ec t s ;^and 

the upper and lower surface flow fields are coupled through a 
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diaphragm region in the plane of the wing. An iterative 
solution is used to arrive at the non-uniform flow solution 
for both non-lifting and lifting cases. Final results for a 
subcritical case, a swept tapered wing at M«»=0.8, show that 
the method converges in three iterations and gives excellent 
agreement with experiment at o^=0° and 2°. Recommendations are 
made for development of a procedure for routine application. 


DEVELOPMENT OF THE METHOD 


Two approaches were investigated for the development of 
a method for coupling linear theory lifting surface solutions 
together to solve the non-linear thick wing problem. In these 
approaches, called the diaphragm and edge singularity methods, 
the upper and lower surfaces of the wing were represented by 
separate non- interfering lifting surfaces whose induced external 
flow fields were forced to match off of the wing but only in 
the same plane of the wing. The boundary conditions required 
that the upper and lower pressure coefficients and the normal - 
wash be equal in the interface region. The differences in the 
two approaches was embodied only in the means by which the 
external flow field conditions were met. The inadequacy of 
either method led to the development of a combined approach, 
the modified diaphragm method, which maintained the best 
features of both. 

In the following paragraphs, the details of the method 
development will be discussed along with results of the 
numerical experimentation. 


Diaphragm Method 


In this method, lifting surfaces on which the free surface 
boundary conditions are met, are added in such a manner as to 
surround the wing being considered. The example shown in 
Figure 1 is the configuration used in this study. The leading 
and trailing edge diaphragms had three and two chordwise 
pressure functions, respectively, as compared with eight on 
the wing. Four spanwise functions were used on all surfaces. 
Although a tip diaphragm should have been included, it was 
eliminated for sake of economy. The solution process was 

identical^ to thab~Tbf^~three wing^~ wi th ~mu tual :in terf erencev^ 

The computer program used was that developed under Contract 

NASl-12399 (ref. 1). 
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For non-lifting cases, only the upper or lower surface 
need to be considered, liie boundary condition met is zero 
downwash on the diaphragm and surface thickness slopes on 
the wing. For lifting cases, both the upper and lower sur- 
faces must be considered. The conditions met in the diaphragms 
are that the downwash and pressure for both upper and lower 
surfaces are equal in these regions. Flow tangency is satis- 
fied on the upper and lower surfaces of the wing for the 
corresponding thickness slopes plus or minus the angle of 
attack, twist, or camber. 

The equation applicable to the non-lifting case is a 
simple interference aerodynaxalc matrix equation 


j^AwP I AwwJ J ^2, J 


( 1 ) 


where 


* A»w 


Awd I Aww 


ao , 
ax 


= aerodynamic sub matrices for downwash 
induced on the diaphragm by the 
diaphragm and wing pressures 

= aerodynamic sub matrices for downwash 
induced on the wing by the diaphragm 
and wing pressures 

= pressure function coefficients for 
the diaphragm and wing 

= surface thickness slope distribution. 


For lifting cases, both upper and lower surface equations 
are written as 



with the conditions 


(3a) 

[w."| = -[w.*-] (3b) 

From equation 3a, a simplification is possible since 

[a:.] = [ Ai.] 

Thus from equations 2, 3b, and 4, the lifting equations become 



The terms and are the surface thickness coordinates as 
modified by angle of attack, camber, or twist. Although 
equation 5 can be reduced by the mamber of rows and colxamns 
corresponding to , it is more convenient to retain this 

form for purposes of iteration. If core storage becomes 
critical, then the reduced form can be used as follows: 



For purposes of this study, equation 5 was used throughout. 
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The effect of the diaphragm is illustrated in Figure 2 
where a non-lifting case is considered for the wing configura- 
tion shown in Figure 1. The airfoil section is a NACA 64A010 
profile normal to the 45° swept quarter chord axis of the wing. 
Experimental pressure distributions on both upper and lower 
surfaces are available in ref. 2 for a range of Mach numbers 
and angles of attack. The solid line is the midspan solution 
with leading and trailing edge diaphragms for the surface 
thickness slopes. The flow is uniform over all surfaces at 
Mjt = Moo= 0.8. The dashed line is the same solution with the 
diaphragms removed. The symbols are experimental data. In 
comparing these results with experiment, it is clear that the 
diaphragm concept provides a significant improvement, particu- 
larly toward the trailing edge. The leading edge region, 
however, is not greatly improved due to the fact that very 
high order functions are still needed to represent the surface 
thickness slopes near the leading edge where they become 
singular for a blunt nose airfoil. 


Edge Singularity Method 

The two- fold purpose of this method is to achieve the same 
result as the diaphragm method with less cost and storage re- 
quirement, and improve the leading edge representation. The 
concept is basically the same with exception that the diaphragm 
regions are replaced with singularities along the edge of the 
wing. The boundary conditions are: 

(1) Stagnation flow at the leading edge (which is 
satisfied by a line source along the leading edge); 

(2) A downwash discontinuity at the leading edge 
between the upwash and mean camber line slope 
(which is satisfied by a line doublet also 
along the leading edge) ; 

(3) Flow tangency along the upper and lower surfaces 
(as satisfied by conventional pressure distri- 
butions) ; 

(4) Zero pressure coefficient differential along 
the trailing edge and strearawise tip; and 

(5) Zero d ownwash^if f e r en tial^j iis t ^ft^ bf^he 
trailing edge and outboard of the strearawise tip. 
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O □ Measured lower and upper 

surface pressures (ref, 2) 

Uniform flow solution with diaphragms 

as shown in fig. 1 

Uniform flow solution, wing only 


Figure— 2, -—Iff ect -o f-^the_diaphragms_ on_ the uniform 

flow solution for surface thickness slopes 



Items 4 and 5 are met with either doublets along the trailing 
edge (one on the upper and one on the lower surfaces) or one 
additional pressure term in both the upper and lower surface 
pressure functions. 

An illustrative example for application of the leading 
edge line source is shown in Figure 3. The line source is 
placed at an appropriate point from which stagnation flow is 
achieved at the leading edge. The strength of the line source 
is determined for the free stream conditions and planform 
geometry, and is held constant throughout the analysis for 
all values of of . The body surface generated by the line 
source, 2 ,, as indicated by the dashed line, is subtracted 
from the airfoil coordinates, 2u and £i,. This step yields 
the modified upper and lower surface coordinates, and e 1, 
which no longer have singular slopes at the leading edge, 

A more specific example is shown in Figure 4 where the 
airfoil section (NACA 64A010 profile swept at 45°, ref. 2) 
is treated with the leading edge line source. The modified 
surface, ^ j , now appears more like a sharp edged 

airfoil. Likewise, the modified slopes no longer tend toward 
a strong singularity. The value of the leading edge radius 
parameter, /?tg= 0.005c, was chosen from graphical considera- 
tions . 

Equations for applying the leading edge line source model 
are developed and presented in the appendix. These equations 
were coupled with the diaphragm method to determine the effect 
of varying the parameter, A diaphragm solution was 

developed at M**=0,8 for the wing shown in Figure 1 where the 
Mach number distribution for non-uniform flow was obtained from 
the experimental data in ref. 2. The effect of varying in 
this solution is shown in Figure 5. The solution for the 
graphically determined value, /?*,£• = 0.005c, is the best of 
those shown, however, - 0,004c would have been slightly 
better. More important, these results imply that the leading 
edge line source coupled with thin wing lifting surface solu- 
tions appears to provide a realistic tool for solving 3-D thick 
wing problems. 
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Figure 3.- Decomposition of a blunt nose airfoil section 







A study was conducted on the application of leading and 
trailing edge doublets to a non-lifting case in order to 
determine the feasibility of the edge singularity method. 

The problem solved was identical to that shown in Figure 5 for 
Rix = 0.005c (solid line), with the exception that leading and 
trailing edge doublets were used in place of the diaphragms 
shown in Figure 1. Two configurations were considered as shown 
in Figure 6. The second configuration differed from the first 
only in the placing of the control points for satisfying the 
trailing edge Kutta condition as shown in the figure. Since 
the upper and lower surface flow fields were S3mimetric, zero 
downwash was the only condition that had to be satisfied at the 
leading and trailing edge external control points. Figure 7 
shows the results for configurations 1 and 2 as compared with 
the diaphragm solution from Figure 5. The most interesting 
effect is the drastic change of the leading edge pressures as 
a result of changing only the trailing edge control point 
location. The use of the leading edge doublet seems to degrade 
the improvement gained with the leading edge line source. 
Finally, the most obvious change is at the trailing edge where 
configuration 1 yields a solution that is similar to that 
obtained without the diaphragm (dashed line. Figure 2). 
Configuration 2 over corrects the problem and further degrades 
the leading edge pressures. 

The use of line doublets to replace the diaphragms has the 
advantage of simplicity, however, it does not appear to be 
adequate. Since the doublets are being used to account for 
flow discontinuities at the edges, it is obvious that the 
solution cannot be rigorous. Logarithmic singular type pressure 
functions would provide the correct characteristics in the edge 
regions but they would greatly complicate the method. 


Modified Diaphragm Method 

The diaphragm method discussed above had the disadvantage 
of not being able to easily treat regions near the leading edge 
on blunt nose airfoils. The edge singularity method had the 
disadvantage of requiring the use of sophisticated singularities 
rather than simple line doublets to meet the edge boundary con- 
ditions. The combination of the leading edge line source with 
the diaphragm method which was used to demonstrate the effect 
of varying (Figure 5), however, appears to be the logical 
ap^oach^ Thi^”mbdifi^~di^hra'g'm”method" maximizes^ the advan- — 
tages and minimizes the disadvantages of both. Although the 
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for non- lifting cases 
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R t.K = 0,005c 

)laa = 0.80 

oL = 0 . 0 ° 

O □ Measured lower and upper 

surface pressures (ref, 2) 

Nominal diaphragm (fig. 1) 

Singularity configuration 1 (fig. 6) 

Singularity configuration 2 (fig, 6) 


Figure 7,- Comparison of edge singularity and diaphragm 
solutions— with— the-experimental- -M^ -distributions 




diaphragm still exists in this method, the total number of 
additional unknown quantities as compared to the edge singular- 
ity method is nominally two in the chordwise direction and one 
in the spanwise direction (if a tip diaphragm is used). This 
is a small price to pay for keeping the. basic simplicity of 
the diaphragm method. 

The equation for the non-lifting modified diaphragm method 
is the same as equation 1 with exception of the right hand side 


A 06 1 A®''*' 


J 

“ 1 

Aw7 y 'a'wT 
_ — 

- 

1 - 4 


1 

1 

/ 


where 

MTw s (&s -r ©; ) (g 

The angle ©5 is the angle of the surface given in the 
appendix 

o, = 

The angle ©V is the flow angle induced by the residual wash, 
normal to the ^ surface and is defined as 


where 


V 

— « f z. “** 0 *1 

L z + cy'-OM; J 


( 9 ) 


It should be noted that ©v is based on the total, local 
tangential velocity, , rather than the value U which is 
induced only by the leading edge line source and the free 
s tream flow ( see eq uations A20) . 
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The equation for the lifting modified diaphragm method 
is likewise similar to equation 5 with exception of the right 
hand side: 


2 /} O0 1 Apv4 j Aowi 


r-'-'i 


® 1 

A wp ' A ! o 

' 

»» ^ 1 
aw 

• r 

"T 1 

i O 1 Aytft 


♦ ^ • 




where 







•f" 



( 10 ) 


(11a) 


c *- 


(11b) 


Although is the same on the upper and lower surfaces, 
will be different as a result of the different velocity fields. 

The next section will be devoted to verification of the 
modified diaphragm method as an economical means for obtaining 
converged solutions for finite thickness wings in the high 
subsonic flow regime. 


VERIFICATION OF THE MODIFIED 
DIAPHRAGM METHOD 


Studies performed with the final method will be sijmmarized 
in this section for both non-lifting and lifting cases in sub- 
critical flow. Convergence of the method will be demonstrated 
and the results compared with a solution obtained with 
distributions taken from experimental data. The effect of 
adding a simple boundary layer model will be shown. Also, 
the importance of accurately determining the surface slopes 
as well as using adequate trailing edge diaphragms will be 
demonstrated. In all solutions obtained, iPtr = 0.005c was 
used. 
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Convergence of the Solution 

Shown in Figure 8 are three iterations for a solution 
to the non-lifting problem which has been the subject of the 
previous section. The case was for oC =0° and Moo =0.80 
which the diaphragm used is shown in Figure 1. The first 
iteration was obtained with a uniform distribution on both 
the wing and diaphragms, equal to M*o=0.80. The distribu- 
tion calculated from the uniform flow solution was used as the 
input for the second iteration. The third iteration was 
obtained with a distribution produced by the second solu- 
tion. The second and third iterations shown in Figure 8 
clearly demonstrate that the solution has converged. 

Figure 9 shows a comparison between the third iteration 
solution and that obtained with the M/ distribution determined 
from experimental data. The differences in used, as indi- 
cated by the theoretical and experimental Cp distributions, 
are much greater than the calculated values of Cp. Hence, 
it is concluded that something may be lacking in either the 
slope distributions used or the diaphragm representation. 

Convergence is shown in Figure 10 for the same wing in 
the same flow with OC=2°, The first iteration was obtained 
with the Mj£ distribution calculated on the third iteration 
above for oC=0°. The second and third iterations in Figure 10 
were obtained in the same manner previously described. In this 
case, the differences noted above are more emphasized in the 
comparison of converged theory and experiment. 

An attempt was made to demonstrate convergence for the 
same wing at Mo© =0.96, ^=0° where the flow became locally 
supersonic. Problems were encountered that produced erratic 
results for iterations which had supersonic flow. The nature 
of the results indicated that a possible error existed in the 
steady flow portion of the computer program used in this study. 
The same program, however, was successfully used to solve 
steady flow problems with mixed subsonic and supersonic flow 
and imbedded shocks (ref. 1). These cases were for rectangular 
wings, which did not require use of certain parts of the pro- 
gram, The suspected parts were examined for possible errors, 
however, nothing was found that could cause what was observed. 
Thus for purposes of this study, the mixed flow case was 
abandoned , 


19 



0 . 

0.005c 

0.80 

0 . 0 ° 

Measured lower and upper 
surface pressures (ref. 2) 

First iteration (M^ = = uniform) 

Second iteration from first 

iteration) 

Third iteration (M^ from second 

iteration) 


Con ve rg enc e’HD f ~the ~theo r e t ica h's o lu t ion~f oroL~=0 
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Possible Sources of Error 


Several effects were examined as possible sources of error 
in the theoretical solution. Because of the closeness in agree- 
ment between solutions obtained with experimental values of 
and theoretical values, the experimental distribution was used 
in these studies. In doing so, the calculated solution should 
agree with experiment if all of the errors were removed. 

Boundary layer effects .- Shown in Figure 11 is the effect 
of including a simple ramp type boundary layer displacement 
thickness, The slope increment was determined from the 

value of i * at the trailing edge obtained from the turbulent 
boundary layer formula (ref. 3) 


The slope was assimed constant at 'bS 0.022, hence, 

the results in Figure 11 only indicate order of magnitude 
and general characteristics of the boundary layer effect. 

This model also will not account for local separation. 

Since the effects of ^ * as shown in this illustration do 
not enhance the agreement between calculated and measured Cp's, 
it is felt that this is not the major problem. 

Effect of errors in measuring surface thickness slopes . - 
The surface thickness slopes were measured again more accurately 
and some differences were found as indicated in Table 1. 


TABLE 1. SURFACE THICKNESS SLOPES 

FOR AN NACA 64A010 AIRFOIL 
AT 45° SWEEP (REF. 2) 


X 

c Original Remeasured 


0.0338 

0.1305 

0.2772 

0.4539 

0.6369 

0.8013 

0v9251 

0.9915 


0.100 

0.046 

0.022 

-0.007 

-0.032 

-0.042 

-0^043 

-0.044 


0.100 

0.046 

0.019 

-0.005 

-0.032 

-0.042 

-0.043 

-0.043 
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^ T.F. = 0.005c 
= 0.80 
06 = 0 . 0 ° 

O □ Measured lower and upper 

surface pressures (ref. 2) 

Surface thickness slopes only 

; Surface thickness slopes plus 

estimated X* slopes 


-1-hi — ff e ut-o f — e s tima t ed— b ounda ry— lay er^d is pi ac emen t 
thickness, S*y on solutions obtained with 
experimental distribution. 



Although the differences were small, the resulting calculated 
Cp's shown in Figure 12 indicate that their impact is signifi- 
cant. With exception of near the trailing edge, the calculated 
Cp's now show a remarkable agreement with measured values. 

Thus, it is concluded that correct measurement of surface 
thickness slopes is a very critical necessity. 

Effect of the trailing edge diaphragm . - As was discussed 
in the evaluation of the edge singularity method, how the 
trailing edge is treated has considerable influence in the 
total solution. Figure 13 shows the configuration with a 
doubled trailing edge diaphragm with three chordwise functions. 
The new diaphragm had the same chord and aspect ratio as well 
as number of pressure functions as the leading edge diaphragm. 
The result of using the larger diaphragm is shown in Figure 14. 
Although the change was not large, the solution was improved 
in the trailing edge region. It is concluded that the use of 
diaphragms needs further study in order to determine the 
relative importance of the parameters involved. 


Final Predictions for oC=0° and oC=2° 

On the basis of what was learned in the above studies, 
final solutions for and oC=2° were assembled with: 

Mjj distributions obtained from the third iteration; enlarged 
trailing edge diaphragm; and remeasured surface thickness 
slopes. 

The results at and oC=0° are shown in Figure 15 

for three span stations, A comparison is made between the 
modified diaphragm method (final solution) , uniform flow over 
a single surface, and experiment. As can be seen, the new 
method is far superior to the uniform flow solution with a 
single surface. The computer costs on a CDC 6600 are shown 
in the following table. 

CP sec 


Single Wing, Uniform Flow 13.5 

1st Iteration 23.2 

2nd Iteration 23.2 

3rd Iteration 23,2 

4th Iteration (expanded diaphragm) 24,5 

Iteratioh^^lutioh^o^al ~94T 1 
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'n = 0 

RlE = 0.005c 
Moo = 0*80 
oC = 0.0° 

O D Measur( 


Measured lower and upper 
surface pressures (ref. 2) 

Original surface thickness slopes 
Remeasured surface thickness slopes 


Figure 12 . - Effec^t^f “the~ac“curacy^of~the“measured— surface 
thickness slopes on solutions obtained with 
experimental M» distributions. 
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Figure 13 .- Geometry and downwash point array for a 

diaphragm solution for the wing of ref. 2 
with enlarged trailing edge diaphragm. 
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Figure 15.- Final non-linear transonic flow solution 
for the wing of ref. 2 at 0C= 0° and Mb*- 




Although the total cost of the iterated solution is not large, 
it could be reduced significantly by automating the procedure. 
Also, if the third iteration were made with the expanded dia- 
phragm, the fourth iteration would not have been necessary. 

The over-prediction at the inboard span station in 
Figure 15 is felt to be due to the effect of wind tunnel wall 
boundary layer. Since the model was a semi-span model mounted 
on the wall without any boundary layer removal, the effect 
could be significant. It was stated in ref, 2 that the esti- 
mated displacement thickness was about 0,5 in. which should 
have a negligible effect on a model with a semi-span of 46.7 in. 
The pressure isobars shown in ref. 2 at oC=2°, however, did not 
verify this assumption since they had a tendency to become 
parallel to the wall near the leading edge. 

The mid-span station shows excellent agreement with 
experiment since it was the least influenced by possible wall 
boundary layer effects and absence of a tip diaphragm. The 
outboard station also shows excellent agreement which was 
apparently due to the fact that the leading edge diaphragm tip 
vortex was in opposition to that shed by the wing. Thus, it 
was felt that large tip diaphragms may not be necessary, 

A comparison is made in Figure 16 for oC=2° between mid- 
span pressures calculated with the original and the remeasured 
thickness slopes. Both solutions were obtained with the 
distribution from the third iteration and the large trailing 
edge diaphragm. The effect of the correction at ot =2° is about 
the same as that seen at oC =0° in Figure 12. 

Figure 17 shows the final results at three span stations 
for Mctt=0.8 and ot =2°. The same comparison is made again to 
show the superiority of the modified diaphragm method. The 
differences between the theories are even more striking at 
06=2°. The computer costs on a CDC 6600 are given in the 
following table. y 

CP sec 


Single Wing, Uniform Flow 13.5 

1st Iteration 41.3 

2nd Iteration 41.3 

3rd Iteration 41.3 

4th Iteration~(expanded diaphragm) —4373. 

Iteration Solution Total 167.0 
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As with the o6=0° case, if the third iteration were made with 
the expanded diaphragm the fourth would not be needed. 

The discussion on disagreement between theory and experi- 
ment given for OC=0° is applicable to ^=2°, It is interesting 
to note that the agreement improved slightly for the mid-span 
station whereas it was degraded forward and improved aft along 
the inboard station. The tip station results showed a degrada- 
tion on the upper surface and an improvement on the lower sur- 
face. Hence, the tip diaphragm may become more important with 
increasing angle of attack or lift. 


CONSIDERATIONS FOR UNSTEADY FLOW CALCULATIONS 

An obvious question that arises when considering the 
calculation of unsteady flow perturbations on the upper and 
lower surfaces is whether or not the solutions can be con- 
sidered as independent in the case of non-sjmunetric flow 
fields. Or, in terms of the modified diaphragm method, would 
it be necessary to obtain a complete upper and lower surface 
solution with diaphragms for unsteady perturbations or would 
simple uncoupled single wing solutions for each side suffice. 

The answer to this question cannot be realized without obtain- 
ing both solutions to the same problem and comparing the results. 
However, some insight can be gained on the basis of what has 
already been done. 

A comparison of the results for ct=0° and oC=2° in 
Figures 15 and 17 shows that the trailing edge pressure coef- 
ficient changes very little at each span station for a change 
in eC. Likewise, the pressure distribution in the diaphragm 
changes very little. Thus, for small changes in lift, the 
presence of the diaphragm is of second order importance in the 
case of subcritical flow at low angles of attack. 

Magnus and Yoshihara (ref. 4) have calculated non-linear 
unsteady pressures on a two dimensional airfoil oscillating 
in mixed transonic flow with an embedded shock. For a wide 
range of reduced frequencies, k=0 to k=5, their results also 
show that the trailing edge pressure coefficient remains un- 
changed at about the steady state value throughout the cycle 
of oscillation. The flows are vastly different on the upper 
and lo we r su rfac es wh ich gives further justification for 
assuming that the flows may be tr ea t ed ' ihd^ehd eh t ly7 
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Thus, it is felt that good results can be obtained with 
decoupled solutions as long as the correct upper and lower 
surfaces Mach number distributions are used. Coupled solutions 
are needed, however, to determine what increase in accuracy is 
gained with what increase in cost. Also, since larger matrices 
will be required which are complex, computer core storage avail- 
ability becomes a limiting factor for the coupled solutions. 


CONCLUSIONS AND RECOMMENDATIONS 


A study has been conducted in order to determine the 
feasibility of calculating non-linear steady flow solutions 
at transonic speeds with the non-uniform flow linear theory 
aerodynamic method of ref. 1. Two approaches were studied 
which included a diaphragm method and an edge singularity 
method. The objective of both approaches was to couple 
together the linear theory aerod 3 mamic solutions for the 
upper and lower surfaces of a wing with finite thickness. 

The investigation of the two methods led to the development 
of the modified diaphragm method which was a combination of 
the two original approaches. 

The modified diaphragm method was demonstrated to be a 
powerful tool for obtaining high subsonic aerodynamic solutions 
for wings with blunt nose airfoils, A leading edge line source 
was used to remove the stagnation flow characteristics from the 
solution which left only a thin wing problem to solve on the 
upper and lower surfaces. The non-uniform aerod)niamic method 
of ref, 1 was then used to solve the thin wing problems where 
the upper and lower surface solutions were coupled only through 
external diaphragm regions surrounding the wing. Free surface 
boundary conditions were satisfied on the diaphragms. 

The method was found to converge rapidly to very accurate 
solutions when used in an iteration scheme in high subsonic 
flow which was subcritical. The case involved a swept tapered 
wing with a s}mimetric airfoil in M«o=0*8 flow at OC=0° and 2°, 
Starting with a uniform flow solution, a converged non-uniform 
flow solution was obtained in two more iterations. Then start- 
ing with a solution obtained with the converged flow field for 
«C=0°, a converged solution was found for oC=2° in two more 
iterations . 
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The leading edge line source and the trailing edge dia- 
phragm had the greatest influence on the aerod 5 mamic solutions 
in the leading and trailing edge regions, respectively. The 
accuracy of the surface thickness slopes was found to be a 
very critical item in the total solution. The result of in- 
cluding an estimated boundary layer displacement thickness was 
inconclusive but the effect was found to be of second order 
importance. 

In an attempt to demonstrate convergence on the same wing 
at Mco= 0*96, OC=0°, some erratic results were encountered in 
iterations which had supersonic flow. The nature of the re- 
sults indicated that an error might exist in the steady flow 
portion of the aerod 5 mamic computer program. The program was 
examined for possible errors, however, nothing was found that 
could cause what was observed. Thus, for purposes of this 
study, the mixed flow case was abandoned. 

. The effect of coupled upper and lower surface solutions 
for unsteady flow was considered. It was felt that good 
results can be obtained with decoupled solutions as long as 
the correct upper and lower surface Mach number distributions 
were used. Coupled solutions are needed, however, in order to 
access the differences that are introduced, if any, by neglect 
ing the effect. 

In conclusion it is recommended that the aerod 3 mamic 
computer program of ref, 1 be modified so as to have the 
capability to calculate non-linear transonic steady mean 
solutions by the use of the modified diaphragm method. The 
following tasks are recommended as a means of achieving this 
capability; 

1, Debug the section of the program responsible 
for the error encountered. (The problem is 
felt to be associated with coordinate trans- 
formations in the steady part of the program 
which are made for swept surfaces in super- 
sonic flow. Hence, the difficulty was not 
noted for rectangular wings.) 
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2. Modify the program to use the more correct 
form of the downwash-pressure function 
integral equation as suggested in ref. 1 

where Uxy and are the velocities at the 
downwash and integration points, respectively, 
and K( ) is the kernel function. 

3. Develop a practical means for achieving 
highly accurate surface thickness slopes. 

. (As an example, surface thickness and camber 
coordinates could be input at 20 points along 
each downwash chord and interpolated with 
either a local cubic or one dimensional spline 
to obtain the slopes.) 

4. Verify the modified diaphragm method for 
mixed flow over both swept and rectangular, 
wings . 

5. Develop an automated iteration procedure 
which is terminated by either a convergence 
criteria or input limits on number of iterations. 

6. Investigate the effect of a diaphragm on 
unsteady aerod 3 mamics. 

7. Compare unsteady predictions with experimental 
results for wings oscillating about both lifting 
and non-lifting steady flow conditions. 

8. Develop an executive capability to assemble 
the generalized forces for upper and lower 
surface solutions, 

9. Conduct a seminar on the capability and 
utilization of the final computer program 
at NASA Langley Research Center. 
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The resulting computer program should provide the capability 
to calculate steady mean flow solutions, oscillatory pressure 
distributions, and generalized forces for wings in mixed 
transonic flow. The effects of surface thickness, twist and 
camber distribution, mean angle of attack, aerodynamic inter- 
ference, and aeroelastic deformation can all be accounted for. 
The flow fields may be mixed subsonic and supersonic flows 
with or without embedded shocks. It appears that the capa- 
bility of the program will be primarily limited only by the 
imagination and ingenuity of the user. 

General D 3 mamics Corporation 
P. 0. Box 748 

Fort Worth, Texas 76101, February 29, 1976 
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APPENDIX 


LEADING EDGE LINE SOURCE 


Calculation of the Velocity Components 

Referring to Figure Al, let the strength per unit length, 
), of the line source satisfy the stagnation condition 
at the leading edge for the velocity component LUcos^i.;. 

The flow potential at point x,y,a due to an infinitesimal 
line segment at /,■»?,$ is 


where 





R 


Cx-j)* + /a^ r** 


The total flow potential,^, is thus 






dl 


(Al) 


(A2) 


where integration is performed along the source line over the 
length Lg. 

Referring to Figure A2, the following geometry trans- 
formations are made to the source line coordinate system: 


X 

- 

X, 

+■ 

X* Cos ■*’ 

y' y/N 

(A3a) 

y 

= 

y. 

+ 

y* cos JtuEK. - 


(A3b) 



X, 

+ 

COS + 

TJ' S lt^~AmL£lt. 

(A3c) 


s 

Xi 

+ 

■>7 'cos -A- It IS. - 

S s/f^ 

(A3d) 




Figure Al.- Leading edge line source configuration 



Figure A2,- Coordinate systems for a 
leading “edge^line' ~^u e 



Since 

()C-j) = (x-x, - 

Cv - ) s ( y. y, - 17* cosJ^lbk) 

and 

= “><• “ +> 5 *f y-y, - v'cos-/i 4 «ie)* 

•f /«'*'C^ “5)*“ 

= Cx-x,y + /^^£cy-y,)‘ + C*-5)‘J 

- 2 j^C X -x«) 5/jw jVux /S^( y~Y, ) QOS A-uaJ '*l * 

+ CO s J '7'* 


Let 

y4 s S/iy*-4.t«<t + /S^cos^JliSii 

B = -aCfx -*,) s/mM-u« + f ‘fy-yOcos-'i-ittJ 
C = Cx-O*- + ^>-[’Cv-y,V * Ca-a,)*J 

Then 

/?«• = A r‘ ^ c 


(A4a) 

(A4b) 

(A4c) 

(A5) 
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Letting 7}t(y) be parabolic. 


>r\(iri*) s </o «4 '7" ♦ </» y* 


equation A2 becomes 


'^*5' LAv'^ * ffV* * 


or 


s 4- d, T, dx Iz 

For Aio, 


(A6) 


(A7) 


2 - ' -gj g *z«/g^ I 
- V5T •^1 Z,«y,' B *zsfA% I > ^ 


(A8a) 


■ ^ h-c^) - ’-rrn^f)] . - ‘ " «»> 


I 0 I B I 


i S'- » 44c 


(A8c) 
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^ - \/3r - 1 r.] (A8d) 

r’-dr,' ^ 

f f'i'- If) ^ - (y.'-iD/jET 

^ AC 

4A 

where 

X , = fk; + c 

s A y»* + BVt ■► c 

If Xi-^O, >/^ 5 O and if Xa<0 ,/X» 5 0. 

For AaO, 



T, - 
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The line source velocity components are calculated with 
a finite difference step as follows: 


(AlOa) 
(AlOb) 
(AlOc) 

where /^K s AV s ^3: “ is the average value of the 

leading edge radius for the line source being considered. 

The coefficients, <it and are determined by satisfy- 
ing the leading edge stagnation condition at the three locations 

y/ 5 y/ + 0.067 ^y*-y.O 

% « y/ +o.5-ooCy/“y/) 

y,' = y/ + 0.9 53 (yi-y/) 








(9,0 

Cos Juieit 

fZu cyi') 

Cos 

CoS 


U - ^ £^(x» 4X, v,a)_ 4fx,y, + U., 

= la = ^C^Cx.y.a+Az)- 4C».y,a)] 


The velocity component for which stagnation is achieved is 

etc - U«* cos-A-tfix c:os -^4« 

The Mach number used in computing ^ with equation A7 and A8 
or A9 is . 
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Calculation of the Surface Zs 

The 2j surface used in this study is defined as that 
surface produced by a two dimensional line source of constant 
strength as required to produce a stagnation point at , 

Thus as RtA changes along the span of the wing, the surface 
likewise changes. 

The appropriate equation relating and X as shown in 
Figure A3 is 


= TT i W(|.) J 


(All) 


for which no direct solution exists for s . 

Approximations have been derived which are applicable in the 
three regions shown in Figure A3, i.e., as , jxl a , 

and the remainder. 

For the first region as X , the inverse tangent 
term may be approximated as 

* -<* _ ( ffu 

^ 3 X* 


IT 


Substitution into equation All yields 




- - 4 ^] 


(A12) 


for ~~ X ~ 
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For|X|->o , equation All is rewritten as 

z, = X ir - |i] 


(A13) 


At Xso , » TTi^ijt/z 9 thus the tangent term may be 

approximated as 



z. 


Now 



which yields 

- TT /ire X 

— /r "ES (ai4) 

for I xl ^ ^ - 

The remaining regions, <X < ~o»i Kls and 

X >o. i fitc are treated with a combined approximation and 
iteration scheme. Starting with equation A13, the following 
approximation is made for as 

?s 

JrRue 


TTPi.6 


= / - e 


- T» C < ♦ /?t^ } 
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Figure A3.- Equation summary for calculating Z 



The coefficient yi is chosen so that 

^ X S - J?LC , O > 
At X = o , 2* r TTRut/Z , thus 


J- s 
2 “ 


I- £ 




Now equation A13 becomes 


^ “ kijt 


0>^93l 

Rl6 


2,^ i Xtcui^Tl- )J 


or 


2s. = X 


t<w% IT 


-0.4931 ( I fe) 


J 


(A15) 


The i^^ approximation is calculated as 

2<. s W Fu r f - 4r tc^~' 1 

*- J (A16) 

which is iterated until | 2sc - / <. i?i.«/iooo . 

Slopes are calculated in all regions except Ixl ^ o,lFu as 


f _ 2^6 — Zs 

Wx “ x»-xT 


where 






wFut 

/ooo 


O . 9 8 F 


and 


X4 






(A17a) 


(A17b) 


(A17c) 
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If 2^4 >.987T^4fi , let 

H s TT /?U5 £ » - 

thus 


2s 


TT Rue 




for which the slope is 



(A18) 


for Zsi 2 .99 . In the region lx| ^ 0.| , the slope 

is calculated as 



Rc6 


4 / 1 *" 

' 4 


Jt. 


which completes the 


equations for slopes of Rs • 


(A19) 


Velocity Components on Zg 

Since the aerod 3 mamic solution is performed along the 
surface Ss , the normal and tangential velocity components are 
needed. Referring to Figure A4, the velocity components in 
the X and ar directions are u. and w* and the tangential and 
normal velocities are Si and^. 
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Figure A4.- Tangential and normal velocities 
on surface Zs 

I 


The total velocity is 

where U. and UT are defined in equations (AlOa) and (AlOb) , 
and the angle is 


(A20a) 


©V - ‘ ) 


(A20b) 
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Thus W, and 5r are 


C3L s V Cos ( ©V - 
^ = V 5ii^ Cfiv - &s^ 

where is obtained as 


(A20c) 

(A20d) 


<^s 


( 


i£/) 


from equations A17, A18, or A19. 


(A20e) 


Generally, 5?“ is non zero since 2^ is defined for two 
dimensional flow. Also, due to the use of the local linear!- 
zation concept, the velocities are calculated according to the 
local Mach niamber, which will differ from those calculated 
with 


Exact Is en tropic Relations JJsed 
in Computing Cp, M^ and 

For a given value of U. as calculated from the above 
equations for the leading edge line source, the pressure 
coefficient, Cp , is calculated as 



(A21) 


The local Mach niimber, , is calculated for a given Cp as 



_j r ^ f 

r-i Lo ♦ 



(A22) 


The tangential velocity, U^, for a given is calculated as 


n . MgU^ r 1 

Moo - 


(A23) 
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